Gluon Propagators in Maximally Abelian Gauge in SU(3) Lattice QCD 
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In SU(3) lattice QCD, we study diagonal and off-diagonal gluon propagators in the maximally 
Abelian (MA) gauge with U(1)3XU(1)8 Landau gauge fixing. These propagators are studied both 
in the coordinate space and in the momentum space. The Monte Carlo simulation is performed on 
16'' at /J=6.0 and 32* at /3=5.8 and 6.0 at the quenched level. In the four-dimensional Euclidean 
space-time, the effective mass of diagonal gluons is estimated as Mdiag — 0.3GeV and that of off- 
diagonal gluons as MoS — IGeV in the region of r = 0.4 — l.Ofm. In the momentum space, the 
effective mass of diagonal gluons is estimated as Mdiag — 0.3GeV and that of off-diagonal gluons as 
Moff — IGeV in the region of p < l.lGeV. The off-diagonal gluon propagator is relatively suppressed 
in the infrared region and seems to be finite at zero momentum, while the diagonal gluon propagator 
is enhanced. Furthermore, we also study the functional form of these propagators in momentum 
space. These propagators are well fitted by Z/ {p^ -\- m?Y with fit parameters, Z,m and v in the 
region of p < 3.0GeV. From the fit results and lattice calculations, all of the spectral functions of 
diagonal and off-diagonal gluons would have negative regions. 

PACS numbers: 12. 38. Aw, 12.38.Gc, 14.70. Dj 



I. INTRODUCTION 



Quantum chromodynamics (QCD) is the fundamental 
gauge theory of the strong interaction based on quarks 
and gluons. There are a variety of nonperturbative 
phenomena in low energy QCD such as color confine- 
ment and chiral symmetry breaking. These nonperturba- 
tive phenomena have been much studied [lH3: however, 
there are still open problems in QCD. In particular, it 
is thought that the confinement mechanism, how quarks 
and gluons are confined, is a considerably difficult prob- 
lem. 

Nowadays, there are several confinement scenarios in 
various gauges. In the Landau gauge, confinement is re- 
lated to the (deep) infrared behavior of the gluon prop- 
agator and the ghost propagator from the scenarios sug- 
gested by Kugo and Ojima, Gribov and Zwanziger [^-ly]- 
These propagators have been much studied from the ana- 
lytical frameworks such as Schwingcr-Dyson equations 
and functional renormalization group equations @ and 
lattice QCD calculations [9l-fl]|. 

On the other hand, the maximally abelian (MA) gauge 
has mainly been investigated from the viewpoint of the 
dual-superconductor picture [l^ - l26| . which is one of 
confinement scenarios suggested by Nambu, 't Hooft 
and Mandelstam [23|. This picture is based on the 
electromagnetic duality and the analogy with the one- 
dimensional squeezing of the magnetic flux in the type-II 
superconductor. In this picture, there occurs color mag- 
netic monopole condensation, and then the color-electric 
flux between the quark and the antiquark is squeezed as a 
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one-dimensional tube due to the dual Higgs mechanism. 
From the viewpoint of the dual-superconductor picture 
in QCD, however, there are two assumptions of Abelian 
dominance [28l . [29| and monopole condensation. Here, 
Abelian dominance means that only the diagonal gluon 
component plays the dominant role for the nonperturba- 
tive QCD phenomena like confinement. 

The various lattice QCD Monte Carlo simulations show 
that the MA gauge fixing seems to support these assump- 
tions [T5l - [25j . In fact, the diagonal gluons seem to be 
significant to the infrared QCD physics, which is called 
"infrared Abelian dominance". Thus the difference be- 
tween the diagonal gluon propagator and the off-diagonal 
gluon propagator in the infrared region seems to be sig- 
nificant. 

From this viewpoint, the gluon propagators in the MA 
gauge h ave been investigated in SU(2) lattice simulations 
lillii,!!^] and in the SU(3) lattice simulation [H. From 
these lattice simulations the off-diagonal gluons do not 
propagate in the infrared region due to the large effec- 
tive mass Mod ~ IGeV, while the diagonal gluon widely 
propagates. In addition, the study of Schwinger-Dyson 
equations also supports the infrared Abelian dominance 
from the behavior of the scaling solutions, i.e., as mo- 
mentum goes to zero, the diagonal gluon propagator is 
divergent and the off-diagonal gluon propagator is van- 
ishing 

The aim of this paper is to investigate the gluon prop- 
agators in the MA gauge of SU(3) lattice gauge the- 
ory. In Secini we briefly summarize the definition of 
MA gauge with U(l)3xU(l)8 and some properties of lat- 
tice gluon propagators in this gauge. In Sec lIIIi we show 
the method of estimating the effective mass in coordi- 
nate space and in momentum space. We study the gluon 
propagators in coordinate space in Sec lIVi and their func- 
tional form in momentum space in Sec|Vl Section IVII is 
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devoted to the summary. 



II. SU(3) FORMALISM AND GLUON 
PROPAGATORS IN MA GAUGE WITH 
U(l)3xU(l)8 LANDAU GAUGE 

Using SU(3) lattice QCD, we calculate the gluon prop- 
agators in the MA gauge with the U(l)3xU(l)8 Landau 
gauge fixing. In the MA gauge, SU(3) gauge symmetry 
is partially fixed, and only U(l)3xU(l)8 Abelian gauge 
symmetry remains. Accordingly, while diagonal gluons 



Al{x){a 



behave as Abelian gauge fields, ofi'- 



diagonal gluons A°(x)(a ^ 3,8) behave as U(l)3xU(l) 
charged matter fields. In the MA gauge, to investigate 
the gluon propagators, we use the glu on fields extracted 
directly from the link- variables (23ll3l|. Here, we analyt- 
ically extract gluon fields A^j^(x) from the link- variables 
U^{x) = e^o-gA^ix) ^]^}^ ^j^g lattice spacing a and the 
gauge coupling g as follows: 



A^{x) = —^l{x)l.nUi{x)n^{2 



lag 



(1) 



where Ln is the natural logarithm defined on com- 
plex numbers, U^{x) the diagonalized unitary matrix of 
U^{x), and ^^{x) the diagonalization unitary matrix. 



U^ix) = %{x)U^{x)^l{x). 



(2) 



We can obtain U'^{x) and Vl^{x) analytically by solving 
a cubic equation with real coefficients using Cardano's 
method j23l |. 

The MA gauge fixing is performed by the maximization 

of 

4 

where H = (r3,T8) is the Cartan generator. In contin- 
uum limit, agA^{x) — > 0, this gauge fixing is given by 
the minimization of 



^ I d''xAl{x)Al{x). 
05^3.8' 



(4) 



Thus, this gauge fixing corresponds to minimizing the 
off-diagonal gluons under the gauge transformation. 

In this gauge fixing, there remains U(l)3xU(l)8 gauge 
symmetry. In order to study the gluon propagators, we 
fix the residual gauge. After the Cartan decomposition 
for the SU(3) link- variables as C/p(x) = M^(x)u^(x) 
with u^{x) ^ ^^{sH'^)T^+s\x)T-) eu(l)3xU(l)8 and 
M^{x) = e*5:„^3,8«°(-)r° eSU(3)/U(l)3xU(l)8, the 
residual gauge fixing is performed by the maximization 
of 



^u(i)L = X! X! tr[u^(a;)]. 



(5) 



In our calculation, this procedure corresponds to "mini- 
mal" MA gauge, and a random Gribov copy is taken in 
the Gribov region [2^ . After gauge fixing completely, we 
study the diagonal (Abelian) gluon propagator as 

G^';s(^ - y) ^ ^ E {Al{^)K{y)) , (6) 

a=3,8 

and the off-diagonal gluon propagator as 

G°^{x-y)^ \ J2 {A;ix)A:{y)) . (7) 

a#3,8 

Note that G'^^^{x—y) and G°^(x— y) are expressed as the 
function of the four-dimensional Euclidean distance r = 
\J {x^ — j/p)^. When we consider the renormalization, 
these propagators are multiplied by an x-independent 
constant, according to the renormalized gluon fields ob- 
tained by multiplying a constant renormalization factor. 

Furthermore, on the periodic lattice of Li x L2 x L3 x 
Li, we consider the momentum-space gluon propagator, 
which is defined as 



a=3,8 



(8) 



(9) 



where p and p are defined as 



27r??.„ 2 . f Pua^ 



(10) 



with 



0, 1, 2, . . . , — 1, and Af^ {p) is defined as 



^m(^5) = E^"''''"'""''"^m(^)- 



(11) 



Note that when the local Landau gauge fixing is satisfied 
as 9pA°(x) = 0, this condition is expressed by 



p^i"(p) = (Landau) 



(12) 



in momentum space. In continuum limit, agAfj^{x) — > 
0, the global gauge fixing condition to maximize i?u(i)L 
in Eq. ([S]) reduces to the local U(l)3xU(l)8 Landau 
condition df,A1^{x) =0 (a = 3,8). Therefore, Eq. ^ 
is satisfied for the diagonal components in momentum 
space. Then, the diagonal gluon propagator has only the 
transverse component T'^^^s{p^): 



Gt^Hp) 



T^'^s(p^). 



(13) 



On the other hand, the off-diagonal gluon propagator 
consists of two components corresponding to longitudinal 
and transverse components. In the continuum limit, the 
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local gau ge fi xing condition for the off-diagonal gluons is 
given by |2l| 



[H, [D^,[D^,H]]]=0 



(14) 



with the covariant derivative Z)^, and this condition is 
different from the Landau gauge fixing ((T^)- Therefore, 
the off-diagonal gluon propagator has the longitudinal 
component even in the continuum limit. Then, the off- 
diagonal gluon propagator is expressed by 



p. 



S,.~^] T°^{p') + (15) 



p. 



with the longitudinal component L°^{p'^) and the trans- 
verse component T°^{p^). These two functions are de- 
rived from the off-diagonal gluon propagator as follows: 



L°^ip') 



,2 GfM. 



P- 

3 



ip) 



PfiPu ^off 



p^ 



G%{P) 



(16) 
(17) 



In this way, we also study the transverse function of the 
diagonal propagator, and the transverse and longitudinal 
functions of the off-diagonal propagator in momentum 
space. 

The Monte Carlo simulation is performed with the 
standard plaquette action on the 16^ lattice with /3 =6.0 
and on the 32^ lattice with /3 =5.8 and 6.0 at the 
quenched level. The lattice spacings a are determined 
so as to reproduce the string tension a ~ 0.89GcV/fm. 
At /3=5.8 and 6.0, the lattice spacings a are estimated as 
a ~ 0.186fm and 0.152fm, respectively |32|. All measure- 
ments are done every 500 sweeps after a thermalization 
of 10,000 sweeps using the pseudo heat-bath algorithm. 
We mainly use 50 configurations for the 16^ lattice and 
20 configurations for the 32^ lattice at each /?. The error 
is estimated with the jackknifc analysis. 



III. PROPAGATOR IN THE PROCA 
FORMALISM 

Next, we investigate the effective gluon mass. We start 
from the Lagrangian of the free massive vector field 
with the mass M ^ in the Proca formalism, 

C = i(a^A,-c),A^)^ + iA/2^^^^, (18) 

in the Euclidean metric. The scalar combination of the 
propagator G^^{r; M) can be expressed with the modi- 
fied Bessel function Ki{z) as 

G^^{r;M) = {A^{x)A^(y)) 



Jk-{x-y) 



(2^)4 

(2^' 
3 M 



Jk-{x-y) _ 



k^ + NP NP 
1 



4,2 ,-Ki{Mr) + —j\x-y). (19) 



In the infrared region with large Mr, Eq. 
to 



3VM e 



-Mr 



2(27r)i ri ' 



(1191) reduces 



(20) 



using the asymptotic expansion, 
Ki{z) : 



^ -.V- r(|+n) 



^^nird-n) (2^)"' 



(21) 



for large Re z. Here T is the gamma function. Then, from 
the slope of the lattice QCD data of In (r^/'^G^^{r)) , the 
effective gluon mass is estimated. 

On the other hand, the propagator in momentum space 
Gf^^{p;M) can be expressed by 



Gm^(p;M) 



1 



Pp.Pu 

AP M2 ■ 



(22) 



Therefore, the longitudinal component L{p^) and the 
transverse component T{p^) in the Proca- formalism 
propagator are regarded as 



Hp') - np') - pi ' 



hp' 



(23) 



Thus, by comparing the lattice QCD data of L{p'^) and 
T{p'^) with Eq. (|23p . the effective gluon mass can be 
estimated. 



IV. ANALYSIS OF GLUON PROPAGATORS IN 
MA GAUGE IN COORDINATE SPACE 

We show in FigH] the lattice QCD result for the di- 
agonal gluon propagator G^^'''(r) and the off-diagonal 
gluon propagator G°^(r) in the MA gauge with the 
11(1)3x11(1)8 Landau gauge fixing on 32^ lattice with /3 — 
5.8 and 6.0. At the long distance, the off-diagonal gluon 
propagator is largely reduced, while the diagonal-gluon 
propagator takes a large value. The infrared abelian 
dominance is thus found in the MA gauge. 

To evaluate the infrared abelian dominance quantita- 
tively, we estimate the diagonal and off-diagonal gluon 
masses in coordinate space from the slope on the loga- 
rithmic plot of r^/^f^diag^^-) g^j^fj ^3/2(^off^j.^^ respectively. 

In Figlll we show the logarithmic plot of f'^^^G°^(r) and 
^3/2(^diag|^^^ as the function of the distance r in the MA 
gauge with the U(l)3xU(l)8 Landau gauge fixing. Note 
that the gluon-field rcnormalization does not affect the 
gluon mass estimate, since it gives only an overall con- 
stant factor for the propagator. We summarize in Table I 
the effective off-diagonal gluon mass Mes and the diago- 
nal gluon mass Mdiag obtained from the slope analysis in 
the range of r = 0.4-0.8fm at /J =5.8 and 6.0 on 32^. The 
off-diagonal gluons seem to have a large effective mass of 
A/off — l.lGcV, while the diagonal gluon mass seems to 




Gjj; ,p=6A32* . 
Gj'«,p=5.8,32* . 
g;'' .p=5.8.32' . 



FIG. 1: The SU(3) lattice QCD results of gluon propagators 
G5l)r^(r) and G'^{r) (top), and their logarithmic plots (bot- 
tom) as the function of r = {x^ — y^)^ in the MA gauge 
with the 11(1)3 xU(l)8 Landau gauge fixing in the physical 
unit. The Monte Carlo simulation is performed on the 32* 
lattice with /3 = 5.8 and 6.0. 



be Afdiag — 0.3GeV. This means that only the diagonal 
gluons , in the MA gauge propagate over the long 
distance and the infrared abelian dominance is found in 
the MA gauge. This result approximately coincides with 
our previous work [23j . 



TABLE I: Summary table of conditions and results in SU(3) 
lattice QCD. The off-diagonal gluon mass Mod and the diag- 
onal gluon mass Mdiag are estimated from the slope analysis 
of r^/^G°^(r) and r^^'^Gf;^'^{r) for r = 0.4 - l.Ofm at each 13. 
In the MA gauge, the off-diagonal gluon mass is estimated as 
Mofi ~ l.lGeV, while the diagonal gluon mass estimated as 
Mdia, =i 0.3GeV. 



lattice size 


P 


a[fm] 


A/off [GeV] A/d 


ag[GeV] 


32* 


5.8 


0.152 


1.1 


0.3 




6.0 


0.104 


1.1 


0.3 



Here, these effective masses are to be considered to give 
the exponential damping of the correlation, and they are 
not simple pole masses. Actually, as will be shown in 



U.l 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



r[fm] 

FIG. 2: The logarithmic plot of r^/^G°^,(r) and r^^'^Gf;^'^{r) 
as the function of the distance r in the MA gauge with the 
U(l)3xU(l)g Landau gauge fixing, using the SU(3) lattice 
QCD with 32* at /3 = 5.8 and 6.0. 



ScclVK. the functional form of the gluon propagator does 
not indicate a simple massive propagator with a definite 
pole mass. 

Note that, because of Mos ^ Afdiag, the diagonal- 
gluon propagation is dominant at the long distance as 



GfJr) 



(r) 



gdiag 



g-(A/oft-A/d.ag)r ^ ^24) 



and, accordingly, the diagonal-gluon contribution be- 
comes decisively dominant in the infrared region. In this 
way, the quantitative difference between the diagonal and 
the off-diagonal effective mass leads to the qualitative 
difference for the long-distance physics. This situation is 
similar to the pion contribution to the nuclear force. At 
the long distance, the nuclear force is well described by 
only one-pion exchange, and the contribution from other 
heavy mesons can be neglected. In fact, the long-distance 
physics is dominated by light fields. 

Next, we also consider the functional form of diago- 
nal and off-diagonal gluon propagators. In our previous 
work [i^l, the off-diagonal gluon propagator is well de- 
scribed by the four-dimensional Euclidean Yukawa func- 
tion as e^"""/r with the mass parameter m in the region 
of r = 0.1 — 0.8fm. In FigI31 we show the logarithmic 
plot of rG°^{r) and rGf}l^^{r) as the function of the dis- 
tance r in the MA gauge with the U(l)3xU(l)8 Landau 
gauge fixing. The logarithmic plot of rG°^(r) seems to 
be linear. From the fit result by the four-dimensional 
Yukawa function in the region of r = 0.1 — l.Ofm, the 
mass parameter m lies in l.SGeV for off-diagonal gluons. 

On the other hand, from FigOl the diagonal gluon 
propagator is not well described by the four-dimensional 
Yukawa function at least in short distance. In Sec|Vl we 
consider the functional form of the diagonal gluon prop- 
agator in momentum space. 

We comment on the two types of gluon-ficld definitions 
on lattice, and their coincidence on the gluon propaga- 



5 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



r[fm] 

FIG. 3: The logarithmic plot of rG°^(r) and rGf;^<^{r) as the 
function of the four-dimensional Euclidean distance r in the 
MA gauge with the U(l)3xU(l)8 Landau gauge fixing with 
32* at 13=5.8 and 6.0. The solid line denotes the best-fit four- 
dimensional Yukawa function. For rG'^{r), the approximate 
linear correlation is found for Yukawa-fit mass parameter m ~ 
1.3GeV. 



From the infrared behavior, these parts would be finite 
at zero momentum, 

L°^{p') ~ T°^{p^) ^ 0. (25) 

On the other hand, in high momentum, the longitudinal 
part is larger reduced than the transverse one. These be- 
havior indicates that the Proca propagator is not suitable 
for the off-diagonal propagator as the functional form in 
the whole momentum region. 

The diagonal part is larger enhanced in the infrared 
region than the off-diagonal parts and thus dominates 
the infrared nonperturbative phenomena. This suggests 
the infrared Abelian dominance. At zero momentum, 
we cannot determine from the infrared behavior whether 
the diagonal part is finite or not. As momentum goes 
to zero, the diagonal propagator seems to be still finite, 
however, for the definite argument in the deep infrared 
region, more careful analysis with a larger lattice size 
would be needed. 



tor. In this paper, the gluon fields ^^(a;) are extracted 



from Ufj.{x) 



On the other hand, with the 



Cartan decomposition of U^{x) = Mfj^{x)u^{x), one can 
take another definition of the diagonal gluons extracted 



'■"■9 E„=3 



! -^^ and the off-diagonal glu- 



In the 



from (x) = e 

ons extracted from M^{x) = e^^s i^a^s.s 
continuum limit, these two definitions of gluon fields co- 
incide. In our lattice calculation, even with these def- 
initions of gluons, we obtain almost the same results 
for both diagonal and off-diagonal propagators. In Ap- 
pendix 1^1 we also estimate the diagonal gluon mass in 
the Cartan-decomposition formalism. 



V. ANALYSIS OF GLUON PROPAGATORS IN 
MA GAUGE IN MOMENTUM SPACE 

A. Gluon propagators in MA gauge in momentum 
space 

In this section, we investigate the momentum-space 
propagators in the MA gauge in SU(3) lattice QCD. In 
Fig|4l we show the transverse and longitudinal compo- 
nents of off-diagonal propagator and the transverse com- 
ponent of diagonal propagator. Here, we calculate them 
with the two different lattice size, 16"' and 32"' at /3 = 6.0. 
As for the diagonal gluon propagator, we only adopt 
32"', because its behavior largely depends on the lattice 
size [2^, while off-diagonal gluons have small volume- 
dependence. Note that, in MA gauge with U(1)3X U(l)8 
Landau gauge, there is almost no the longitudinal diag- 
onal component, which is checked in Appendix IB] 

In low momentum {p ^ IGeV), the behavior of the 
longitudinal part is similar to that of the transverse part. 



''II, 




p [GeV] 




p [GeV] 

FIG. 4; The logarithmic plot of {p'^),L°^ (p^) and 
ydiag|-p2-j ^ ^j^g function of the momentum p = {p^p^y^^ 
with 16* (top) and 32* (bottom) at /?=6.0. The diagonal 
gluon propagator is calculated only on 32* due to large vol- 
ume dependence. In infrared region, the diagonal gluon prop- 
agator is largely enhanced, while the longitudinal and trans- 
verse parts in the off-diagonal gluon propagator are relatively 
suppressed. The infrared Abelian dominance is found also in 
momentum space. 
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B. Estimation of gluon mass in MA gauge in 
momentum space 

First, we estimate the effective mass of gluons by com- 
paring the obtained gluon propagators with the massive 
vector-boson propagator in momentum space. The lon- 
gitudinal and transverse components in the Proca for- 
malism have the common factor, according to a constant 
renormalization factor Z of the renormalized gluon fields, 

^b^) = TiP^) - (26) 

As for the off-diagonal propagator, the transverse com- 
ponent is fitted by p-i^^^fi with the fit parameters, Z and 
M, in the region of p < Pmax = l.lGeV. Fixing the con- 
stant Z ^ we also fit the longitudinal component by pi^j^p 
with a fit parameter M in the region of p < Pmax- 

The diagonal transverse component is also fitted in a 
similar manner to the off-diagonal transverse one. 

We summarize in Table II the fit results in the region of 
P < Pmax = 1.1 GcV. The off-diagonal gluon masses from 
the transverse and the longitudinal parts almost coincide 
as 

Mjff ~ - IGeV. (27) 

On the other hand, the diagonal gluon seems to behave 
with the effective mass MT-„ ~ O.SGeV from the trans- 
verse part result. These results coincide with the coordi- 
nate results in Scc lIVI 

TABLE II: Summary table of the ghion mass estimation in 
momentum space. The off-diagonal gluon mass Moff and the 
diagonal gluon mass Mdiag are estimated from the fit analysis 
by the Proca propagator for p < l.lGeV. In the MA gauge, 
the off-diagonal gluon mass is estimated as Moff — IGeV, 
while the diagonal gluon mass estimated as Mdiag — 0.3GeV. 
These results coincide with the coordinate results. 



lattice size 


Pmax[GeV] A/[GeV] 


Z 






1.1 0.93(2) 


2.7(1) 


0.9 


32* 


1.1 1.00(2) 


2.7(fix) 


2.5 


ydiag 324 


1.1 0.26(2) 


22.0(8) 


1.7 



off-diagonal gluon propagator has not be described by 
the propagator of the Proca formalism, but by the four- 
dimensional Euclidean Yukawa function e~™''/r with a 
mass parameter m. The Fourier transformation of this 
Yukawa function is expressed by 

/ d^xe'P-- 1 ^ , „ . (28) 

J r (p^ -|- m'^y^i'^ 

Thus, we fit the transverse and longitudinal part of 
the off-diagonal propagator and the transverse part of 
the diagonal propagator by 



{p^ -\- to2)3/2 

with the fit parameters Z and m in the region of p < 
Pmax = S.OGcV. Note that Z cannot be regarded as a 
renormalization constant, and, in general, the fit param- 
eters Z in the transverse part and the longitudinal part 
would differ. 

We show the fit results in Figl5]and the parameters in 
Table ITTll As for the off-diagonal gluon propagator, both 
transverse part and longitudinal part are well described 
by the four-dimensional Euclidean Yukawa function, al- 
though, these parameters are a bit different. This result 
coincides with our previous work [2^ and this means that 
the spectral function of off-diagonal gluons has a nega- 
tive region [2^ [s^]. On the other hand, the diagonal 
gluon propagator is not well described with the four- 
dimensional Euclidean Yukawa function. 

TABLE III: The fit results of the off-diagonal parts and 
the diagonal part by the four-dimensional Euclidean Yukawa 
function for p < 3.0GeV. The transverse and longitudinal 
parts in the off-diagonal gluon propagator are well described 
by the four-dimensional Yukawa function with /N ~ 1.4, 
while the diagonal gluon propagator is not well described from 
the large xV^ (~ 2.4). 





lattice size 


Pmax[GeV] 


m[GeV] 


Z 






32* 


3.0 


1.27(0) 


6.08(4) 


1.4 




32* 


3.0 


0.95(1) 


3.22(3) 


1.4 




32* 


3.0 


0.36(1) 


17.7(1) 


2.4 



C. Analysis of functional forms in MA gauge in 
momentum space 

As the fit region increases, gluon propagators cannot 
be described well with the massive vector boson prop- 
agator, pi^j^j-i ■ In fact, in the region of p < S.OGeV, 

X^/iV of the off-diagonal transverse part is larger than 3, 
and that of the diagonal transverse part is larger than 
15. This result seems to be consistent with our previous 
work [23|. In coordinate space, the functional form of the 



More generally, we fit the propagators by more general 
Ansatz, 



(30) 



(p2 + 



with the fit parameters Z, m and v in the region of 
P < Pmax = 3.0GeV. We call this "i/- Ansatz" . We show 
the fit results in Fig 15] and the parameters in Table HVl 
The transverse and longitudinal parts of the off-diagonal 
propagator have different functional forms: The trans- 
verse part has v ~ 1.3—1.4 and the mass parameter 
m ~ IGeV and the longitudinal part has ;/ ~ 1.8 and 



7 



O 1 




FIG. 5: The logarithmic plot of T°^{p^), L°^{p^) and 
ydiag^p2^ as the function of the momentum p = (p^p^Y^^ 
with 32* at /3=6.0. The transverse and longitudinal parts in 
the off-diagonal gluon propagator are well described with four- 
dimensional Yukawa function, while the diagonal part is not 
well described in the infrared region. All parts of gluon prop- 
agators are well described with !/-Ansatz (|30|l in the whole 
region of p < S.OGeV. 



m ~ IGeV. From this result, the transverse part behaves 
Hke a four-dimensional Yukawa function with the mass 
parameter m ~ IGeV. On the other hand, the longitudi- 
nal part docs not behave Hkc it {v ~ 1.8 > 1.5), however, 
the mass parameter takes a similar value, m ~ IGeV. 
Due to this larger the longitudinal part is larger re- 
duced than the transverse one in high momentum. The 
diagonal part has ~ 1.8 and the mass parameter 
0.6GeV. The diagonal part is better described with these 
parameters than the four-dimensional Yukawa function. 
Note that both diagonal and off-diagonal gluon propa- 
gators show 1/ ^ 1, which means that they cannot be 
described by a simple massive propagator with a definite 
pole mass. 



TABLE IV: The fit results of the off-diagonal parts and the 
diagonal part by the //-Ansatz (|30|) for p < 3.0GeV. All these 
parts in gluon propagators are well described by this fit func- 
tion with corresponding parameters. 





lattice size 


Pmax[GeV] 


m[GeV] 


u Z 






32* 


3.0 


1.09(2) 


1.34(2) 4.0(2) 


1.2 




32* 


3.0 


1.21(4) 


1.75(4) 6.0(7) 


1.2 


ydiag 


32* 


3.0 


0.58(1) 


1.76(1) 29(1) 


1.2 



D. Spectral function of gluons in MA gauge 

The spectral function is given by the inverse 

Laplace transformation of the zero-spatial-momentum 
propagator DQ{t). Then we define 



-1 /'C+ioo 
^TTl Jc-ioa 



(31) 



where D^'^it) is defined with the gluon propagator 
G^'''{x) by 



(32) 



Note that D'^^{t) and p^''(w) no longer have Lorentz ten- 
sor structure due to the spatial integral. Corresponding 
to the two components of the gluon propagator. 



Tip') + ^L{p') } (33) 



the zero-spatial-momentum propagator has also two com- 
ponents. 



dpA 



J (2^) 



(2^) 

Note that Drit) is obtained from the space-space correla- 
tor and DlII) is from the time-time correlator. Thus, we 
can define the two components of the spectral function 
as 

p'^'^iiu) = (5^" - S^^S"^) pt{uj) + Sf'^S^^pLiio), (35) 
where 



Pt,l{^) 



2m 



dte'^^DT,L{t). 



(36) 



Pt,l(w) denotes the spectral function of transverse and 
longitudinal gluons. In the same way, we can define the 
corresponding effective mass. 



MT,Lit) = In 



DT.Ljt) 
^T,L(t + l) 



(37) 
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In the Landau gauge, as t increases, the effective mass in- 
creases. If the spectral function p{uj) is non-negative, the 
effective mass M{t) should be a monotonically decreas- 
ing function. In fact, if M{t) has an increasing part, p{uj) 
should have negative region [s^ . 

If the functional form of some part (longitudinal or 
transverse part) of the gluon propagators is well de- 
scribed by the j/-Ansatz, (pa_|,f„:^).^ > 0)j the cor- 
responding zero-spatial- momentum propagator D^{t) is 
obtained by 



> 
o 





Mj(l),offdiiigonal i — * — i 




1 

, i I 



dp4 

(2^)' 

Z 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 
t [fm] 



{p. 



^/^^T {v) \2in 



- 

/-1/2 



K,-i/2imt) (Rei>0),(38) 



where K^_i/2 is the modified Bessel function. Therefore, 
the corresponding effective mass Myit) is expressed by 



M,{t) = In 



= In- 



i/2{'mt) 



(39) 



The modified Bessel function is reduced to Ku-i/2{'m''t) 



(mt) 



-1/2, 



for large mt, and thus 



M,{t) ~ In 



— l^ — m 



. 1 

m — {i> — 1) —. 



(40) 



This indicates that, if the functional form of the propa- 
gator is well described by the i/-Ansatz M^{t) in- 
creases for 1/ > 1 as <: increases. Due to this, the spectral 
function would have a negative region. 

In fact, we calculate the off-diagonal effective masses 
and the diagonal effective mass. There are two off- 
diagonal effective masses corresponding to the longitu- 
dinal and transverse parts, however, there is only the 
transverse part in the diagonal effective mass, due to 
(9pAJ^ ~ (a = 3,8). In FiglBl we show these effective 
masses in SU(3) MA gauge with U(1)3X U(l)8 Landau 
gauge. We calculate the off-diagonal effective masses on 
16"' with P — 6.0 and the diagonal effective mass on 32'* 
with P = 6.0 due to the large dependence of the lat- 
tice size. For the diagonal gluon and the off-diagonal 
transverse/longitudinal gluons, each effective mass M{t) 
is found to be an increasing function of t at least in small 
t region, which leads to existence of negative region of 
the spectral function of each gluon component 33 1. 



VI. 



SUMMARY AND CONCLUDING 
REMARKS 



We have studied the gluon propagators in the MA 
jauge with the U(l)3xU(l)8 Landau gauge fixing using 



> 
o 



ML(t),offdiagonal 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 
t[fm] 

My(t), diagonal ■ — 



> 
o 



0.4 

0.3 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 
t [fm] 

FIG. 6: The effective mass plot of tlie longitudinal and trans- 
verse parts of off-diagonal giuons, Mhit) and MT{t), and the 
diagonal gluon. The off-diagonal effective masses are obtained 
on the 16* lattice with /? = 6.0 using 200 configurations, and 
the diagonal effective mass is obtained on the 32* lattice with 
/3 = 6.0 using 40 configurations. Each effective mass M{t) 
seems to be an increasing function of t. 



the SU(3) lattice QCD both in coordinate space and in 
momentum space. The Monte Carlo simulation is per- 
formed on the 16* lattice with (3 ~ 6.0 and 32* with /3 
=5.8 and 6.0 at the quenched level. 

In coordinate space, we have calculated the Euclidean 
scalar combinations of the propagators G'^p(r) in the di- 
agonal and the off-diagonal gluons and estimated each 
gluon mass from the linear slope of the logarithmic plot 
of r^/^Q^^(^r) We have found that the off-diagonal gluons 
behave as massive vector bosons with the effective mass 
A'/off — 1 GeV for r = 0.4 — 1.0 fm, while the diagonal 
gluons behave as lighter vector bosons with Mdiag — 0.3 
GeV for r — 0.4 — 1.0 fm. Due to the larger off-diagonal 
gluon mass, the off-diagonal gluons cannot mediate the 



interaction over the large distance as r ^ -^off ' while 
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the diagonal gluons can propagate in this region. Such a 
behavior would lead infrared Abelian dominance in the 
MA gauge. 

Furthermore, we have investigated the diagonal gluon 
propagator and ofF-diagonal gluon propagator in momen- 
tum space. To our knowledge, the analysis of the gluon 
propagators with the MA gauge of SU(3) lattice in mo- 
mentum space is the first study. In the MA gauge with 
11(1)3 X U(l)8 Landau gauge, the diagonal propagator 
has the only one component and the off-diagonal prop- 
agator has two components, the transverse part and the 
longitudinal part. In the infrared region, the diagonal 
propagator is enhanced, while the transverse and the lon- 
gitudinal parts in the off-diagonal propagator show the 
similar suppressed behavior. In this way, the infrared 
Abelian dominance is also found in momentum space as 
in coordinate space. Furthermore, at zero momentum, 
the two components of the off-diagonal propagator are 
finite and relatively suppressed and the diagonal prop- 
agator seems to be finite and enhanced. However, as 
for the diagonal propagator in the deep infrared region, 
more careful analysis with a larger lattice size would be 
required for the definite argument. 

We have also estimated these gluon masses from these 
components of diagonal propagator and off-diagonal 
propagator in momentum space by compared with the 
Proca formalism in the region of p < 1.1 GcV as in co- 
ordinate space. As a result, the off-diagonal gluon mass 
seems to be Mos — 1 GeV and the diagonal gluon mass 
seems to be Mdiag — 0.3 GeV. These results are consis- 
tent with the analysis of coordinate space. 

In addition, we have also investigated the functional 
form of the propagator in the MA gauge. The gluon prop- 
agators have been fitted by Z/ {p^ -f m-^) with the pa- 
rameters, Z, m and v, in the region of p < S.OGeV. These 
propagators are well described with the functional form. 
The best fit results show that these mass parameters 
mdiag, m^if and m^g have a relation, m^gjm^g » mdiag 
and all of these exponentiation parameters, v, are larger 
than unity. Thus even the functional form of the diag- 
onal propagator is larger enhanced than that of the off- 
diagonal propagator and would be found infrared Abelian 
dominance. Furthermore, from these fit results, in par- 
ticular > 1, all of the corresponding spectral functions 
would have negative regions as expected in our previous 
work [2^. In fact, we have calculated these zero-spatial- 
momentum propagators, I?o(Oj and the corresponding 
effective mass, M{t) = \nDo(t)/Do{t + 1). We expect 
the negative regions of these spectral functions from the 
behaviors of these effective mass. 
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Appendix A: Comparison of Cartan decomposition 
and ordinary decomposition 

In SU(iVc) lattice QCD, the gluon field Aj. is usually 
defined from the link- variable based on [2| 

U^{x) = e'^s^-^") e SU(7Ve), (Al) 

which we call ordinary decomposition (OD). In our pa- 
per, this gluon field is mainly used for the argument 
of gluon propagators. However, in terms of the partial 
gauge fixing in the MA gauge, one can take other defini- 
tion of the gluon field based on the Cartan decompo- 
sition (CD), 

e SU(7Ve)/U(l)^=~' X U(l)"="\ (A2) 

In the continuum limit, these two definitions of gluon 
fields coincide. Actually, in our lattice calculation, even 
with these two definitions of gluons, we obtain almost the 
same results for both diagonal and off-diagonal propaga- 
tors. 

In this appendix, wc compare the diagonal gluon prop- 
agators with OD and CD, which are defined by 

i ^ {A;ix)Atiy)) (OD), 

a=3,8 

i J2 {A;ix)AUy)) (CD). (A3) 

a=3,8 

In Figl?! we show the logarithmic plot of the scalar 
combination of the diagonal gluon propagators with OD 
and CD in coordinate space. This result is obtained on 
32"* with [3 = 6.0. These propagators are almost the same 
behavior. 

Furthermore, we estimate each diagonal gluon mass 
from the slope on the logarithmic plot of ?'^^^G^'^8(r) 
as in Sec lIVI In FiglHl we show the logarithmic plot of 
j,3/2(^diag^j,^ with OD and CD. These gluon masses are 
obtained from the slope analysis in the range of r = 0.4 — 
O.Sfm. In both cases, diagonal gluon mass seems to be 
about 0.3 GeV. Also for the off-diagonal gluons, the CD 
results are almost the same as the OD results. 
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> 




0.001 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 
r [fm] 



FIG. 7: The SU(3) lattice QCD results of Gf;^<^{r) with Car- 
tan decomposition (CD) and ordinary decomposition (OD) 
in the MA gauge with the U(l)3 xU(l)8 Landau gauge fixing. 
The Monte Carlo simulation is performed on the 32'* lattice 
with P — 6.0. The CD gluon propagator almost coincides 
with the OD gluon propagator. 




0.1 0.2 0.3 0.4 0.5 0,6 0.7 0,8 0.9 1 



1- [fm] 

FIG. 8: The logarithmic plot of r^^'Gf;^<^{r) with Cartan 
decomposition (CD) and ordinary decomposition (OD) in the 
MA gauge with the U(1)3XU(1)8 Landau gauge fixing on the 
32* lattice at /3=6.0. The slopes of these plots in the region 
of r = 0.4 — l.Ofm are almost the same. 



Appendix B: The longitudinal part of the diagonal 
propagator 

The lattice global Landau-gauge fixing is given by max- 
imization of 

R[U] =RcTr"^ Uf,{x) (Bl) 

From this global gauge condition, the local gauge fixing 
is derived as 

AR = R[U^'] - R[U^=^] 

= ReTr^[FC/^(x) + [/^(x - fi)V^ ~ (V = 1)] 

< 0, (B2) 

where V £ SU{3) denotes the arbitrary local gauge trans- 
formation on the site x. By expanding the local gauge 



transformation as V = e^^ ^ ~ 1 + iO'^T"^, the local 
gauge formation is obtained, 

Tr J2 - '^m(^) - (Uf^ix - A) - Ulix - A)) } 

= 0. (B3) 
In the Landau gauge, the gluons are usually defined by 

A^{x) = 1 {[/^(a;) - ulix)} - (trace part) . (B4) 

With this definition, the local gauge fixing is expressed 
by the ordinary Landau gauge fixing, 

Tr T- Y.^ iM^) - - A)} = 0. (B5) 

However, with our definition of gluons Eq. ([T|), the local 
gauge fixing condition is given by 

Tr T« J2^, {sin - sin - A)} = 0, (B6) 

where we use precise relation of smA^(x) = 
^ {J7^(x) — Ul{x)}. In small lattice spacing, this local 
gauge fixing corresponds to the ordinary Landau gauge. 
In general, however, this gauge fixing condition does not 
coincide with the ordinary gauge fixing condition, and 
thus, the gluon propagator with our definition of gluons 
in momentum space would have the longitudinal part as 
well as the transverse part. 

The similar situation arises for diagonal gluons in MA 
gauge with U(1)3X U(l)8 Landau gauge fixing, because 
the gauge fixing for U(1)3X U(l)8 part is given by the 
maximization of 

4 

^u(i)L = ^^Iletr[u^(x)]. (B7) 

X /i — 1 

Thus, the local condition is given by 

{smAl^ix) - smAl^ix - fi)} = (a = 3, 8), (B8) 

where A^{x) ~ (A'^ix) , A^{x)) is defined by Ufj{x) = 
(.i<^9T.a=3,8-^'^MT\ This condition (|B8]) is satisfied ac- 
curately in numerical calculation. In this gauge, even 
with small lattice spacing, it is not so trivial whether the 
local gauge condition (jBSp coincides with the ordinary 
Landau gauge due to the large fluctuation of diagonal 
gluons. Note that there are two different definitions of 
diagonal gluons, A1{x) and A'^{x), however, the behav- 
ior of these gluon propagators would be similar as shown 
in Appendix [Al 

In Figini we show the longitudinal part and the trans- 
verse part of the diagonal propagators. The transverse 
part is dominated in whole momentum region, 

2^diag^p2) < T^'^s(p2). (B9) 
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Therefore, there would be the only transverse part in 
the diagonal gluon propagator even with our definition 
of gluons, 

Gf^^p) (^S^, - rd-8(p2). (BIO) 

Note that as for the off-diagonal gluon propagator, 
there are both longitudinal and transverse parts even in 
the continuum limit, 

GfAp) = (-5,. - ^) T°ff(p2) + ^L°ff(p2).(Bll) 



'•iit.^"»».. 




0.1 
0.01 
0.001 

1 2 3 4 5 

p [GeV] 

FIG. 9: The logarithmic plot of T'^'^'^ip'^), L'^'""^{p'^) and 
Gf;^'^{p'^) = 3r'*'="5(p^) + L'*'='*5(p2) as the function of the mo- 
mentum p = (p^p^y^^ in the MA gauge with U(1)3XU(1)8 
Landau gauge fixing with 32* at (3=6.0. In this gauge, 
7;^diagj-p2^ < T'^'"'i^{p'^) is satisfied at least in the region of 
p < S.OGeV 
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